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Strongly correlated electron systems require the development of new theoretical 
schemes in order to describe their unusual and unexpected properties. The usual 
perturbation schemes are inadequate and new concepts must be introduced. In 
our scheme of calculations, the Composite Operator Method, is possible to re- 
cover, through a self-consistent calculation, a series of fundamental symmetries 
by choosing a suitable Hilbert space. 
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The discovery of new materials with a large variety of unusual and unexpected 
properties [ has opened a new era in the physics of Condensed Matter; new theo- 
retical schemes must be developed [ §]. The most important characteristic of these 
new systems is a strong correlation among the electrons that makes inapplicable 
classical schemes based on the band picture. It is necessary to pass from a single- 
electron physics to a many- electron physics, where the dominant part will be the 
correlations among the electrons. Usual perturbation schemes are inadequate and 
new concepts must be introduced. 

Let us consider a certain Hamiltonian 

H = H[<pi(x),...,<p n (x)] (1) 



the set {ifi} denotes band-electron fields. Due to strong correlation effects the prop- 
erties of the original electrons are drastically changed; new excitation modes 
will appear and determine most of the observed properties of the system. It is nat- 
ural to identify a new set of elementary excitations {ipi} as basis for perturbative 
schemes. These excitations, constructed from the original electron set (in this sense 
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we call them composite fields), are created by the interactions among the electrons; 
therefore, their properties will be determined by the dynamics, the symmetries of the 
model, the boundary conditions and must be computed in a self consistent way [ [| . 
This aspect introduces a new richness into the theory and will allow us to realize the 
dynamics in the proper Hilbert space where the physical symmetries are conserved. 
On the other hand, we know from experiments that highly correlated systems ex- 
hibit an incredible variety of behaviors. It would be very hard to describe such a 
complexity using the original fields, unless the exact solution of the model is avail- 
able. The presence of new excitations, composite fields, introduces into the theory 
the possibility to accommodate the multifariousness of experimental properties. 

A theory built on the basis of new excitation modes is, by construction, a self- 
consistent theory, and a procedure must be fixed. In particular, we must answer to 
the following list of questions: 

1. the identification of the fundamental set; 

2. the statistics and the properties of the new fields; 

3. the symmetry and the dynamics in terms of the new fields; 

4. the representation where the new fields are realized. 

We will now try to formulate a scheme of calculations in which an answer to the 
previous questions can be found. Then, by considering a particular model, we will 
present a practical realization of the theoretical scheme. 

The new fields are generated by the interactions among the bare fields; then, 
it is naturally to choose the new set as the one which naturally appears through 
the equations of motion. The evolution of the original fields is described by the 
Heisenberg equation 

d 

i Qi ( Pi( x ) = = Ji [tp(x)] (2) 

Such an equation generates new fields { Ji[<p(a;)]}, constructed as combinations of 
the bare fields. By starting from these fields and by considering the new Heisenberg 
equations 

»j^*(aO = [^(aO,#] (3) 

we generate an infinite hierarchy of composite fields. It is naturally impossible 
to solve the infinite system of equations and some truncation procedure must be 
adopted. Let us consider a n-component field 



0= : 

and let us choose the first n — 1 fields such as 



(4) 



d i+1 
ittMx) = Mx),H] =5>y(-«V)^(z) 1 < j < n - 1 (5) 

01 3=1 
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The n th -field ip n (x) is determined by the field equation of tp n -i(x). The matrix 
7(— iV) is completely determined by the dynamics. Then, we linearize the Heisenberg 
equation by writing 

d 

% ^V X ) = £ ( _?V ) ^( x ) ( 6 ) 

where the eigenvalue or energy matrix e is self-consistently calculated by means of 
the equation 

e(-iV.) i),^(y, *)}) = ({[^(x, 0, #] , V%, *)}) (7) 

The symbol (■ • •) denotes the thermal average. Derivative operators as A(— iV) 
are defined through the relation 

A(-zV)/(x)= /AA(x,y)/(y) (8) 

The rank of the energy matrix is equal to n, the number of components of the 
vector ip{x). When there is translational invariance we can invert Eq. ([^D and it is 
easy to see that 

e ij = {{[^(x),H],ipj(y)}) ({inix^iy)})" 1 = 7 « 

1 < i < n - 1, 1 < j < n (9) 

This approximation corresponds to the n-pole expansion of the Green's function 
where finite life-time contributions are neglected. It has been proved [ || that in this 
approximation the choice (Q) for the composite operators leads to the conservation 
of the spectral moments. In particular, the first 2(n — i + 1) spectral moments for 
the field < i < n — 1] are conserved. This is an important property when we 
recall that the spectral moments are related to the spectral density function of the 
single-particle propagators. Also, as shown in Ref. [[§], the choice (||) leads to an 
equivalence between the n-pole approximation and the spectral density approach 
[ [H, although very different results are obtained when different procedures for the 
self-consistency are used [ |J. 

In general the composite fields will not satisfy canonical anticommutation rela- 
tions and their algebra must be calculated starting from the canonical algebra of the 
electron fields. Owing to this fact, the Wick theorem and the standard perturbation 
schemes cannot be applied. Examples of the new algebra will be presented in the 
second part of this article. 

The properties of the new fields are fixed by a series of parameters which must be 
self-consistently calculated. These parameters are expressed as expectation values of 
composite fields. When the composite fields belong to the set they can be expressed 
in terms of the single-particle Green's function and calculated by a series of coupled 
self-consistent equations. 

However, it may happens that some of the parameters are expressed as expecta- 
tion values of higher-order composite fields that do not belong to the basic set. In this 
case, owing to the approximation considered, the parameters are not strictly bound 
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by the dynamics and there is a freedom in the procedure to fix them. At this level 
the powerfulness of the scheme manifests itself: one can use this freedom to choose 
the right representation. In the construction of a physical theory we must distin- 
guish two levels. On one side we have the microscopic level where we are concerned 
with particles. The basic ingredients are the Heisenberg fields which together with 
the canonical commutation relations describe the dynamics. The physical laws (the 
equations of motion, the conservation laws, the symmetry principles) are expressed 
as relations among the operators. On the other side we have the macroscopic world 
where we are concerned with average values of operators. At the level of observation 
the physical laws manifest themselves as relations among matrix elements, and a 
suitable choice of the Hilbert space must be made. When some approximation is 
introduced the states are not the exact eigenstates of the Hamiltonian; the expecta- 
tion values are also not the exact ones. As a consequence, the relations among the 
operators are generally not conserved when the expectation values are calculated. 
A striking example of this is the violation of the Pauli principle. A convenient way 
to take care of it is to operate in the representation of second quantization where 
the Pauli principle manifests through the algebra. It is known [ [7| that in most 
of the approximation schemes this symmetry is violated when matrix elements are 
considered. Other examples of symmetries will be considered later. 

The point of view adopted in this approach is that we can use the freedom in 
the procedure to fix the self-consistent parameters in such a way to recover the sym- 
metries violated by the approximation. In general, a model exhibits many different 
symmetries and there will be a relation between the number of composite fields 
and the number of symmetries that can be recovered. On the physical ground one 
must choose which symmetries are the most important to be satisfied. In a physics 
dominated by strong electron correlations the Pauli principle plays a crucial role, 
and it is extremely relevant that the related symmetries be treated in a correct way. 
Therefore, in our scheme the attention is firstly put to the Pauli principle; once this 
is accommodated, the attention is devoted to other symmetries. 

As an illustration of the scheme we shall now consider the Hubbard model [ |J. 
In a standard notation this model is described by the following Hamiltonian 

# = E fei -^a) c f (0 c(i) + £/£n t (i) Hii) (10) 

ij i 

c(i), c^(i) are annihilation and creation operators for electrons at site % in the spinor 
notation 




n a = cj.(z) c a (i) is the number operator of electrons with spin a = (t, I) at the 
i th site. jU is the chemical potential and is introduced in order to control the band 
filling n. For a two-dimensional squared lattice and by restricting the analysis to 
first nearest neighbors, the hopping matrix tjj has the form 

iV k 
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a(k) = -[cos(M+cos(V)] 



(13) 



a being the lattice constant. In addition to the band term, the model contains 
an interaction term which approximates the interaction among the electrons. In 
the simplest form of the Hubbard model, the interaction is between electrons of 
opposite spin on the same lattice site; the strength of the interaction is described 
by the parameter U. 

The electron field c(i) satisfies the Heisenberg equation 



d 

i-^c(i) = -fj, c(i) - At c a (i) + U r){i) 



where 



(Tit) 



(14) 



(15) 



is the electron field on the nearest neighbor sites. We see that the dynamics has 
generated the composite field 

r)(i) = c(i) n(i) (16) 

The Heisenberg equation for this field will generate a new higher order composite 
field. The process does not stop and an infinite number of composite fields will be 
generated. By following the procedure mentioned above, we close the hierarchy by 
considering n fields and we construct the vector composite field as described by 
Eqs. (U) and (H). For the specific case, we consider the three-component field 



I Mi) 
M*) 
V Mi) 



(17) 



where 



Mi) 
Mi) 

Mi) 



£(i) = c(z) [l-n(i)] 
rj(i) = c(i) n(i) 
1 
2 



7T U 



~a» ntl (i)c a (i) + c(i) c 1 (i)c(i) 



c\i) c{%) is the charge (/i = 0) and spin (/i = 1, 2, 3) density operator for 



c-electrons. We are using the following notation 



l,a , a" 



•1, a), a being 



the Pauli matrices. The composite fields ( |T7D do not satisfy a canonical algebra. For 
example, for the first two fields 



{r){i),rf(j)} = 



(19) 



5 



F. Mancini and A. Avella 



This field satisfies the Heisenberg equation 







i-Jj(i) = J(i) = (-» + U) V (i) + At tt(z) j (20) 

1 -nir(i) + 4i«(i) -U0(i) + U p(i) 



where 



= ^c(i)c+ a (i)^^(i)-^C«(i)c+(i)«7 M d a[ (i) 

0(z) = -f^ rv(z) c a2 (z) + c a (z) c ta (z) c(i) - c(i) c to2 (z) c(z) + c(z) c fQ (z) c Q (z) 
p(z) = I^^(z)^(z) + c(z)e tQ «c(z) (21) 



According to the method given above, the equation of motion (|2~0"D is linearized 

d 



as 

4^) = £^iMi) ( 22 ) 

where the energy matrix s(i,j) is the 3x3 matrix given by 

e{i,j) = ( { J(0 , ^ (j ) } > E T < { ^ (0 (i) } >^ 1 T (23) 

The subscript E.T. indicates that the anticommutators are evaluated at equal time. 

The physical properties can be described in terms of the thermal retarded Green's 
function 

2 

SUj) = (RUU)^Ui)]) = ^— [ cfkdioe^^-^-^-^ S(k,u) (24) 
x L J/ (27r) d Jn B 

where R is the usual retarded operator and the symbol (• • •) denotes the thermal 
average. By means of the linearized Heisenberg equation (|22|) the Fourier transform 
is given by 

S(k,") = E rffL- (25) 

~ x u - E n (k) + i T] 

where the energy spectra £i(k) are the characteristic values of the matrix e(k), 
determined by the equation 

£ o™(k) = (26) 

m=l 

The characteristic coefficients aj(k) are defined by the following relation 

a n . k (k) = (-) k Tr k [e(k)} < k < 3 (27) 

where Tr^ is the trace of the k order, defined as the sum of the determinants of 
all ( , ) matrices of order k x k which can be formed by intersecting any k rows 
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of e with the same k columns. We note that Tr 3 [e] = Det[e] and the convention 
Tro[e] = 1 is used. The spectral functions are given by 

<r»(k) = -i- £ K l (k) A m (k) (28) 

°nl K J m=0 

where the A n (k) are the 3x3 matrices: 

3 

A"(k) = fl ™( k ) e m ~ n - l {y) I(k) < n < 2 (29) 

m=n+l 

and we put 

b n (k)= n [E n (k) - E m (k)] (30) 

m=l,m^n 

By standard arguments, the correlation functions can be calculated from the 
knowledge of the retarded Green's function. By means of (^) we have 

2 3 

C(i,j) = (m ^(j)) = £ J cPkeiHKi-^-iEnQtM-ti) a n {k) {1 + Tn(k)] 

(31) 

where we put 

r.<k) = t™h(§&) (32) 

We see that the calculation of the Green's function requires the knowledge of the 
normalization matrix 

7(k)=^({^),^ (j )}) ET (33) 

and of the m-matrix 

m(k)=^({j(z),^ (j )}) ET (34) 

where T indicates the Fourier transform. 

These quantities are calculated in appendix and depends on a series of parame- 
ters, that can be so listed: 

1. external parameters as the temperature T and the electron density n = (c^(i) c(i) 

2. model parameters as U and t; 

3. self-consistent parameters that can be calculated in terms of elements of the 
Green's function, as \i and A; 

4. self-consistent parameters expressed as expectation values of composite fields 
out of the basis ([□]), as p, I$ 3 , m° 33 , m%. 
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For the latter a procedure of self-consistency must be fixed. In the Composite Oper- 
ator Method (COM) we take advantage of this freedom and we fix the parameters 
in such a way that the Hilbert space has the right properties to conserve the rela- 
tions among matrix elements imposed by symmetry laws. In a physics dominated 
by a high correlation among the electrons, the first attention should be put to the 
requirement that the approximation does not violate the symmetry required by the 
Pauli principle. Let us consider the correlation matrix fl3T|); when we take equal 
points the algebra leads to the following relations 

(e(<)i7 f (<)) = 

(£(0^(0) = (d°(0 (35) 

(»7(0 **(<)) = -(^(0^(0) 

among matrix elements of the Green's function. These relations constitute a set of 
coupled self-consistent equations which will be satisfied by an appropriate choice of 
the parameters. 

The recovery of the Pauli principle does not exhaust all the degrees of freedom 
and we have place to accommodate other symmetries. An intrinsic symmetry of the 
Hubbard model is the pseudospin SU(2) symmetry, which is nothing also that the 
invariance under the particle-hole transformation. The generators of this transfor- 
mation are given by the total pseudospin operators 

p+ = E(-) ic r« c l« 

i 

P- = £(-yq«c T (.) (36) 

i 

P* = ~£K*)-1] 

i 

These operators satisfy the SU{2) algebra 

[P+, P-} = 2P Z [P ± , P z \ = T P ± (37) 
and the Heisenberg equations 

i|V = ±{2 l i-U)P ± 

= (38) 
Let us consider the thermal retarded Green's function 

P+-(t-t') = (R \P + (t)p-(t') ]) = — / due-***-*) P + ~(uj) (39) 

By means of the equation of motion (j38|) we obtain for the correlation function 

1 (n-Ve-* 2 "-™-*) 



- {P + (t)F-(t')) = x : e ^_ v) (40) 
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This is an exact result which relates the pseudospin correlation function to the 
particle number n and it is a manifestation of the intrinsic symmetry. 

Another important symmetry is given by the conservation of the current density. 
By defining the charge p(i) and current j(z) densities as 

p(i) = ec\i)c(i) (41) 
j(i) = -itea 2 c\i)[V - V]c(i) (42) 

it is immediate to obtain by means of the Heisenberg equation ( |l"4"|) the conservation 
law 

V-j(i) + J^(i) = (43) 



The symmetry content of the algebraic equation P5[ ) manifests at level of ob- 
servation as relations among matrix elements once a choice of the physical space of 
states has been made. Indeed, by defining the causal charge and current functions 

as 



Xab{h3, 

where 



(T[g a (i)g b (i)}) = ^33/ d'kdue^-^—^Xa^Ku) (44) 



p{%) for a = 
9a{i) = { 3x(i) for a = x (45) 
j y (i) for a = y 

we can derive a series of Ward-Takahashi identities connecting current- current, 
charge-current and charge-charge propagators. One of those reads as follows 

iau xoo(k,oj)= [l - e ~ ik * a ] x±o%u) + [l - e~ ik » a ] x y o(K u) (46) 

In the static approximation of the Composite Operator Method the charge, 
current, spin, pseudospin correlation function can be connected to convolutions of 
single-particle propagators. This occurrence is related to a linearized dynamics to- 
gether with the choice of occupation dependent electronic excitations as basic fields 
[ |J. Once these calculations have been performed, Eqs. (|36|), (|4(iD and ( fi"6|) con- 
stitute a set of five coupled self-consistent equations which can be satisfied by an 
appropriate choice of the five parameters p, I^ 3 , m® 3 , m 33 . 

We thus have a scheme of calculations in which it is possible to recover, through a 
self-consistent calculation, a series of fundamental symmetries by choosing a suitable 
Hilbert space. 

Detailed calculations will be presented elsewhere. 
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Appendices 



A. The normalization matrix 



From the definition fl33|) and by means of the canonical algebra for the c-electrons 
it is straightforward to see that for a paramagnetic ground state the normalization 
matrix has the following expression 

J n (k) J 13 (k) 
I(k) = | I 22 (k) J 23 (k) | (A.1) 
/ 13 (k) 7 23 (k) J 33 (k) 



with 



/u(k) = l-£ 

J 13 (k) = A + (p - J 22 ) a(k) 

/ 22 (k) = - (A.2) 

J 23 (k) = -A-pa(k) 

/ 33 (k) = 4 + a(k) i£ 



The quantities introduced in (A.2) are so defined 



e(i)e«)-(^)^ T «; 

A = {t a (i)?(i))-{v a (i)v\i)) (A.3) 



n is the average number of electrons per site; A and p are static correlation function 
between nearest neighbor sites. In particular, the parameter p describes intersite 
charge, spin and pair correlations. In the calculation of I 33 (k) only the nearest 
neighbor contributions have been retained 

({7r(z),7rt(j)})^^/ 3 ° 3 + ^/ 3 c 3 (A.4) 



B. The m-matrix 



At first we note that time translational invariance requires the m-matrix be 
hermitian 

( mn(k) mi 2 (k) m 13 (k) 
m(k) = mi 2 (k) m 22 (k) m 23 (k) | (11/1; 
\ m 13 (k) m 23 (k) m 33 (k) 
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From the definition fl34|) and by making use of the expressions (|20|-f21"|) for the 
source current it is possible to calculate 



(B.2) 



mn(k) 


= -[/i + 4ta(k)]/n-4t/i3 


m 12 (k) 


= -4ta(k)/ 22 -4t/ 23 


mi 3 (k) 


= - + At a(k)] A3 - 4t [a(k) J 23 + I33 


m 22 (k) 


= -(At-C/)/ 22 + 4i/ 23 


m 2 3 (k) 


= -(/x-C/)/ 23 + 4i/ 33 


m 23 (k) 


= "4 + «( k )™33 



In the calculation of m 33 (k) only the nearest neighbor contributions have been re- 
tained 

({js(i),^{j)}) = S ij m° 33 + a ij m^ (B.3) 
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